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A model consisting of oxygen-occupied and -vacant chains is considered, with 
repulsive first and second nearest-neighbor interactions Vi and 1^, respectively. The 
statistical mechanics and the diffraction spectrum of the model is solved exactly and 
analytically with the only assumption Vi >> V2. At temperatures T ^ Vi only a 
broad maximum at (1/2,0,0) is present, while for | 5 — 1/2 |> 1/14 at low enough T, 
the peak splits into two. The simple expression for the diffraction intensity obtained 
for T « V\ represents in a more compact form previous results of Khachaturyan 
and Morris, extends them to all 5 and T/V2 and leads to a good agreement with 
experiment. 



PACS numbers: 74.70. Vy, 61.14.Hg, 64.60.Cn 
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It is now clear that the superconducting transition temperature of YBa2Cu307_5 depends 
not only on 6 but also on the oxygen ordering Therefore, detailed knowledge of this 
ordering in the whole oxygen concentration range is important for an understanding of the 
electronic properties of the system and the superconducting mechanism. 

At low enough temperatures, at least for 1/4 < 6 < 1/2, the oxygen atoms in the basal 
planes are ordered in infinite CuO chains Experimental [|]-|^ and theoretical [f?HTO(| 

evidence favoring alternative types of oxygen ordering is restricted to 6 > 0.6 and 6 < 0.25. 
Concerning the ordering among different CuO chains at T = 0, the experimental situation 
favors structures that are given by the ground state of a one-dimensional (ID) Ising 
model in which the interactions Vn between chains at a distance na, where a is the lattice 



parameter, satisfy the inequality K.+i + K.-i > 2Ki |1C-13|. The low temperature thermody- 



namics of the two-dimensional asymmetric next-nearest-neighbor Ising model (ASYNNNI) 



1^ is also governed by the simplest one- dimensional Ising model {Vi > and = for 



n>l) PJII 



Due to the sluggish oxygen kinetics at low temperatures, the system often does not reach 
a completely ordered state at T=0, as is required by the third principle of thermodynamics. 
In this case, in the range 1/4 < 6 < 1/2, diffuse diffraction peaks are observed 0,0- For 
the largest values of S, and diffraction vectors q = (27r/a)(/i, 0, 0), only one peak in the 
interval 0< h <1, centered at = 1/2 is present. For the smallest values of 6 in the above 
mentioned range, two maxima centered at h = l/2±e are observed. Khachaturyan and 
Morris (KM) explained qualitatively these observations in the range 1/3 < 5 < 1/2, 
assuming a random faulting of the double period ordered structure with 6=1/2 (alternating 
Cu-0 and Cu- vacancy chains 0). This work has been criticized in a Comment [1^ because 



of the restricted composition range of the theory and the assumptions on the interactions 
that would be implicit in the model. In their Reply |T^ KM state that Ref. does not 



assume a particular type of interaction and that, while it would be nice to have a tractable 
analysis for all 6, it was not necessary to establish that random faulting can produce split 



peaks. In another Comment to it was shown that the short-range order implicitly 



assumed by KM, minimizes the free energy of the ID Ising model with repulsive Vi, V2 and 



Vn=0 for 7T, > 2 at T ^ pO|. Vi > 2V2 should hold to stabilize the double period structure 



at 6 = 1/2. While interactions at larger distances are important in determining the ground 



state of the system |T0|-|T3| , it will be shown in this letter that this model is the simplest one 
which leads to a reasonable agreement with the experimentally observed diffraction peaks. 
However, for T the free energy should be minimized by long-range ordered structures 
and this model becomes unrealistic. In addition, the results of KM predict too narrow and 
intense split peaks near 5 = 1/3 when compared with experiment P,[T6[|. Since one expects 
that these peaks should broaden and lose intensity when the temperature is increased, it is 
of interest not only to extend the theory of KM to all compositions, but also to study the 
model at finite temperatures. This task is carried out in the present letter. 

In order to obtain simple analytical results, two cases are considered: a) T >> V2, 
Vi/T arbitrary ||21|, b) Vi » T and any V2/T with V2 < Vi/2. The resulting short-range 
correlations are used as input parameters for the calculation of the diffraction intensity. If 
the system is metastable, these parameters can be thought of as probabilities that are given 
by the preparation method, independent of the free energy of the system, or as equilibrium 
parameters at higher T. While KM obtained two different expressions for the diffraction 



intensity, one for 1/3 < 6 < 4/9 and another one for 4/9 < 6 < 1/2 IT^I, in case (b) 



the elegant method based on generating functions leads to a single expression which 



simplifies those of KM, extends them to all oxygen compositions, and allows for more than 
two consecutive oxygen occupied chains (requiring T 7^ in the ID Ising model). 

Following KM we shall denote by Q ^he Cu-0 chains and by □ the Cu-vacancy chains. 
It is convenient to write the model in the form: 

H = ViJ2riini+i + V2j2^i'^i+2, (1) 

i i 

where = 1 (0) if the ith chain is n(0)- The diffraction intensity is given by [p!^] : 



I{h) = N\U\^6 P{m)exp{t27chm), (2) 



where is the total number of unit cells, fox is the oxygen scattering factor and P{m) is 
the conditional probability that if = 1, also rim = 1- Interchanging Q ^^^^d □ it is easily 
seen that, excluding integer h, I{h) is the same for 6 and 1—6. 

In the following 5 < 1/2 is assumed, 
a) V2 « T 

In this case V2 can be neglected (as will become clearer in case (b)) and the model 
reduces to the simplest ID Ising model. This model describes the low-temperature oxygen 
ordering of the 2D ASYNNNI model [|l3[|l5l and is exactly solvable The quantity that 



determines the free energy and I{h) for each value of 6, is the probability yi of finding a 
pair □□ of two consecutive □ chains: 

yi = 6- 71/2 + [{6 - 7i/2)2 ^ ^2^^^ _ 1/2 ^ 

where 71 = [1 — exp(— Vi/T)]~^. P{m) satisfies the following non-homogeneous difference 
equation: 

P{m) =r]- (3P{m-l), (4) 

where rj = (5 — — 5) is the conditional probability that if n^-i = 0, then = 1- 

Similarly rj — (3 = yi/S = P{1) is the probability that if Um-i = 1 then rim = 1- The solution 
of Eq. (4) with the boundary condition P(0) = 1 is: 

PM=5 + (l-5)(-/3)H, (5) 

and replacing this in Eq. (2) 

/w^A-i/.j- ^'!-/''^;^-''' (6) 

1 + p + 2p coslznh) 

For any 6 and T this expression gives only one peak centered at h = 1/2. If enough 



statistics and Monte Carlo steps are allowed, Monte Carlo results using the ASYNNNI |T5 
at low enough temperatures should converge to this simple expression. Thus, the ASYNNNI 
should be extended to include longer range repulsions in order to explain split diffraction 
peaks |]T0HT3[|, as shown before pO 



In the limit T » Vi |2|, then 71 = T /Vi, yi = 5'^[1 - (1 - 5fVi/T], r] = 5 + 5\l - 
6)Vi/T, (3 = 6{1- 6)Vi/T and: 

I{h) = N\fo, I' 5(1 - 5) [1 - 25(1 - 6){Vi/T) cos(27r/i)] . (7) 

Thus, at high enough T, I{h) is a constant plus a small harmonic term with maximum at 
h = 1/2 and minimum at /i — and /i — > 1 . 

For Vi » T » V2, neglecting exponentially small terms, 71 = 1, = 0, 77 = /5 = 
6/(1-6) and: 

^ N\U\'6il-6)il-26) 
^ ' 1 + 25(1 - 5) cos(27r/i) ' ^ ' 



in agreement with Ref. ||20|. For 5 — > 0, I{h) is small and flat, while for 5 — > 1/2, the second 
member of Eq. (9) gives the extremely narrow peak at h = 1/2 corresponding to the double 
period long-range ordered structure, 
b) Vi » T 

As shown above, for 5 < 1/2 the probability yi of finding a strip □□ is of order 5^(1— 25) 
exp {—Vi/T) and can be neglected if 5 7^ 1/2. This allows to represent any possible structure 
in terms of a sequence of two strips: Q '-'O- The energy per strip is given by V2?/ where 
y is the probability of finding two nearest-neighbors DQ strips. Calling x the probability 
of finding one DQ strip, it is easy to see that the problem can be mapped into the ID 
Ising model with only Vi ^ , already considered in case (a). The mapping is given by the 
correspondence DQ — ^ ^ 6, y ^ yi, V^^Vx- The free energy per chain F is given by: 

(1 + x)F = V2?/ — T[xlnx + (1 — x) ln(l — x) — ylny — 2{x — y) ln(x — y) 

-(1 -2x + ?/) ln(l-2x + ?/)], (9) 

where 1 -|- x is the average number of chains per strip. The average number of □ chains per 
strip is x/{l + x) = 6. Thus: 

x = 5/(l-5), (10) 
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and minimizing F one obtains: 

2/ = x - 7/2 + [(x - 7/2)' + x2(7 - 1)] , (11) 

where 

7 = [1 - expi-V^/T)]-' . (12) 

The probabihty of finding three consecutive Q chains, which is given by 1 + S{y/x — 3) can 
be used instead of T as independent variable. 

P{m) can be determined from P(m — 1) and P(m — 2) from the following reasoning. 
Since P{m) = P{—m), m > will be assumed. The pair of chains m — 2 and m — 1 can 
be in one of the three following states: 1) 00) 2) DO; 3) O'^- The probability of state 
i is denoted pi. Since □□ is excluded, if the third state is realized {p^=l) then P(m)=0. 
In state 1, the chain m — 1 should belong to a O strip and if pi =1, P{m) is given by the 
conditional probability (x — — x) that after a strip O; ^ strip DO follows. Similarly, if 
P2 = 1, P{itl) is the conditional probability y/x that after a strip nO; another of the same 
kind follows. Thus: 

P{m) = pi{x - y)/{l - x) + p2y/x, (13) 
and using P(m — 2) = p2, P{m — 1) = p^, Pi + P2 + Ps = 1, Eq- (13) takes the form: 

P(m) =/?[l -P(m- 1)] -aP(m-2), (14) 

where 

/3 = (x - y)/{l - x) ; a = f3 - y/x. (15) 
Eq. (14) can be solved using the generating function [22]: 



G{z) = J2 PMz"". (16) 

Using Eq. (14) and the boundary conditions P(0) = 1, P(l) = 0, an equation for G{z) is 
obtained, the solution of which reads: 
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r(z) - ^ + (^-^)^ (17) 

By integration in the complex plane it can be shown that 

P{m)=S-Rm+i{zi)-Rn,+i{z2), (18) 

where Rn{zi) is the residue of G{z)/ z'^ at the pole Zj, and Zi and Z2 are the two roots of the 
polynomial az^ + I3z + 1. However, the diffraction intensity can be obtained directly from 
the generating function. Using Eqs. (2) and (16) one has: 

I{h) ^N\fo:c? 8 [G'(exp(i27r/i)) + G'(exp(-i27r/i)) - 1] , (19) 

and after some algebra, this expression is simplified to: 

T(h\ = iV|/o. \^5il-a)il + a-f3) 

^ ' 4acos2(27r/i)+2/3(l + a)cos(27r/i) + (l-a)2 + /52 ^ > 

Eq. (20), together with Eqs. (10) to (12) and (15) describes the scattering intensity for all 
5 < 1/2 and T « Vi_. For 5 > 1/2, 1(h) is given by the same equation with 5 replaced 
by 1 — 5. The condition for the existence of two split maxima is obtained requiring that 
the denominator of Eq. (20) as a function of cos(27r/i) has a minimum in the interval (-1,1). 
One obtains 4q; > /3(1 + a), or equivalently 7 < 7c, where 

7c = 4(1 - 26)/{l - 6) (21) 

Since 7 > 1, split peaks are possible only for 5 < 3/7 — 0.429. For these values of S, the 
simple Eqs. (12) and (21) give the critical temperature above which only one peak exists. In 
the region of compositions and temperatures for which two intensity maxima exist {5 < 3/7 
and 7 < 7c) , their positions are given by a very simple expression: 



hfnax — 2^ arccos 



7(1 - S) 



(22) 



4(l-2(5)_ 

In the limit Vi » T » V2, neglecting terms of order ex.p(—Vi/T) and of order V2/T, 
which do not bring any qualitative change in I{h), one has 7 = T/V2, y = x"^, (5 = x, a = 



and I{h) takes the form of Eq. (8). For T << V2, neglecting exponentially small terms 
one has 7 = 1 and i) for 6 < 1/3, y = 0, a = (3 = 6/{l — 26); ii) for 6 > 1/3, y = 
{36 — 1)/(1 — 6), (3 = 1, a = (1 — 26)/ 6. Case (ii) coincides with the one previously solved 
by KM |T^, since Eq. (14) takes the form of the nonhomogeneous difference equation of 



KM. In both cases, I{h) takes the form (for T « V2 and any 6): 

Uh) = iV|/.. P^(l-2^)|1-3^| 

^ ^ 4(5( 1 -2(5) cos2(27r/i)+ 2(5(1- 5) cos(27r/i) + 1052 + 1-6(5 ^ ^ 

Eq. (10a) (vahd for 4/9 < 5 < 1/2) and also Eq. (10b) of KM (1/3 <6 < 4/9) should reduce 
to the much simpler Eq. (23). Moreover, Eq. (23) extends the results of KM to all values 
of 6, and Eq. (22) with 7 = 1 gives the position of the split maxima. This low temperature 
limit can be described as an uncorrected sequence of i) strips Q □00^^'^^ 1/3 
or ii) strips DQ and □ Q O ^ ^ V^- These strips would become correlated if V3 (and 
V4 for 6 > 1/3) were included in the model. 

For a comparison with experiment, the low temperature limit Eq. (23) is not good enough 
and Eqs. (10) to (12), (15) and (20) should be used. In Fig. 1, the evolution of I{h) with 
temperature is represented. For 6 = 0.364 and T = (7 = 1), I{h) shows two sharp peaks 
as already shown in Fig. 1 of KM . As expected, the peaks broaden and lose intensity. 



keeping the same total area, when the temperature is increased. However, as long as two well 
defined maxima exist, the positions of the peaks do not depend strongly on temperature. For 
7 = 1.2, the result is very similar to one of those obtained by Beyers et al. for 6 = 0.35 p|,p!6|. 
For 6 = 0.25, the experimental peaks are somewhat sharper, suggesting that repulsions at 
larger distances than two lattice parameters are also present |1T0|-[T3|. Other difficulties 
in the comparison between theory and experiment are the possibility of phase separation 
PJT6|] and that a fraction of oxygen atoms always remains disordered [Q, particularly for 
quenched samples. For comparison with the experimental results in quenched samples |T6 



with 6 = 0.27, 6 = 0.35 and 6 = 0.43, 6 is replaced by 55/4 in the theoretical curves, 
following Ref. [^. The corresponding results for 7 = 1.2, shown in Fig. 2, are in good 
agreement with experiment. The intensity for 6 = 0.35 x 5/4 is somewhat higher than the 
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experimental one. 

The agreement with experiment can be improved by adding more interactions. Also a 
quantitative agreement with experiment was obtained postulating that P{2n + 1) = within 
domains, the size of which is adjusted for each 6 [|T6|. However, except in unrealistic limits, 
an analytical treatment of the problem is no longer possible in these cases and a further 
improvement of the present results can only be done at the cost of a loss in simplicity 
and physical transparency. Furthermore, in view of the above mentioned experimental 
uncertainties, the analytical results are already satisfactory. 

I acknowledge gratefully insightful conversations with J. Garces and H.Wio. 
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FIGURES 

FIG. 1. Intensity as a function of diffraction vector for several values of V2/T (see Eq. (12)), 
V\ » T and (a) 5 = 4/11, (b) 5 = 1/4. In (a), the critical temperature above which only one 
peak exists is given by 7 = 7c = 12/7 (see Eqs. (12) and (21)) 

FIG. 2. Intensity as a function of diffraction vector for 7 = 1/2, Vi » T and several values of 5. 
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